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ABSTRACT 

We investigate the detailed processes working in the drift of magnetic fields in 
molecular clouds. To the frictional force, whereby the magnetic force is transmitted to 
neutral molecules, ions contribute more than half only at cloud densities njr <^ 10 4 cm -3 , 
and charged grains contribute more than about 90% at nn ^ 10 6 cm~ 3 . Thus grains 
play a decisive role in the process of magnetic flux loss. Approximating the flux loss 
time ts by a power law ts oc i3~ 7 , where B is the mean field strength in the cloud, we 
find 7 ~ 2, characteristic to ambipolar diffusion, only at nn <^ 10 cm~ 3 where ions and 
smallest grains are pretty well frozen to magnetic fields. At njr > 10 7 cm~ 3 , 7 decreases 
steeply with bh, and finally at bh ~ n^ec ~ a few x 10 11 cm -3 , where magnetic fields 
effectively decouple from the gas, 7 <C 1 is attained, reminiscent of Ohmic dissipation, 
though flux loss occurs about 10 times faster than by pure Ohmic dissipation. Because 
even ions are not very well frozen at njr > 10 7 cm -3 , ions and grains drift slower than 
magnetic fields. This insufficient freezing makes ts more and more insensitive to B 
as 71h increases. Ohmic dissipation is dominant only at nn ^ 1 x 10 12 cm -3 . While 
ions and electrons drift in the direction of magnetic force at all densities, grains of 
opposite charges drift in opposite directions at high densities, where grains are major 
contributors to the frictional force. Although magnetic flux loss occurs significantly 
faster than by Ohmic dissipation even at very high densities as «h ~ ^deo the process 
going on at high densities is quite different from ambipolar diffusion in which particles 
of opposite charges are supposed to drift as one unit. 
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1. Introduction 

Magnetic fields of interstellar molecular clouds are widely believed to have significant effect 
on star formation. For an oblate cloud or cloud core of mass M contracted along field lines to 
some extent, there is a critical value for its magnetic flux <P given by 

$cr = UG 1/2 M, (1) 

where G is the gravitational constant and is a dimensionless constant. A cloud (core) with 
<P < <P cr cannot be kept in hydrostatic equilibrium by the magnetic force alone (this state is widely 
called magnetically supercritical), and one with ^ > ^ cr (magnetically subcritical) can be in 
equilibrium if its expansion is suppressed by external magnetic fields. Applying the virial theorem 
to such clouds, Strittmatter (1966) found between 4.9 and 9.4 depending on the flatness of 
the cloud. From some numerical cloud models Mouschovias & Spitzer (1976) obtained ~ 8.0, 
and Tomisaka, Ikeuchi, & Nakamura (1988) got « 8.3. Li &; Shu (1996) found = 2ir for 
self-similar, singular, isothermal clouds. 

We can define the critical magnetic field strength of the cloud (core) by 

Bcr = ^~ 2 = UG^U, (2) 

where R is the radius of the cloud (core) perpendicular to the mean magnetic field direction 
and U = M/irR 2 is the mean column density of the cloud (core) along field lines. Of course 
B ^ B CT is equivalent to <P ^ ^ cr , where B = <P/irR 2 is the mean field strength in the cloud (core). 
Nakano &: Nakamura (1978) found that isothermal disks penetrated by uniform magnetic fields 
B perpendicular to the disk layers are gravitationally unstable only when B < B cv , where B cr is 
given by equation (|2|) with = 2tt. The critical wavelength, below which the disk is unstable, 
decreases as B / B cr decreases ( Nakano 1988| ). These mean that with perturbations at least the 



disk as a whole can contract perpendicular to field lines when B < B CI . 

While magnetic fluxes of clouds estimated by observations are not much smaller than their 
critical values, or <P/<P cr ~ 1 (e.g., Crutcher 1999] ), magnetic stars with mean surface fields of lkG 



to 30 kG have ratios <&/<P cr ~ 10~ 5 — 10~ 3 ( [Nakano 1983j ). For ordinary stars like the Sun with 
mean surface fields of ~ 1 G, <&/<& cr is as small as 10 -8 . This suggests that cloud cores must lose 
most of their initial magnetic fluxes at some stages of star formation. At what stages and by what 



mechanisms? This is called the magnetic flux problem in star formation (e.g., Nakano 1984 ). 



Ohmic dissipation is too slow to dissipate magnetic fields in molecular clouds of ordinary 
densities. For example, for the ionization fraction 10~ 8 at density 10 5 cm -3 (Figure 1) and the 
length scale of magnetic fields, 0.1 pc, about the Jeans length at this density and temperature 
10 K, we obtain the Ohmic dissipation time 10 15 yr, larger than the age of the universe by orders 
of magnitude. Mestel & Spitzer (1956) found another process of decreasing magnetic flux, which 
is now widely called ambipolar diffusion (sometimes called plasma drift; see Appendix |C] for the 
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terminology) and is much more efficient than Ohmic dissipation in molecular clouds of ordinary 
densities. In this process ions (and electrons), which are well frozen to magnetic fields, drift in 
the sea of neutral molecules together with magnetic fields at a terminal velocity with which the 
magnetic force balances with the frictional force exerted by the neutrals. 

Because some dust grains in clouds are electrically charged and interact with magnetic fields, 
they contribute to controlling the drift of magnetic fields as well as ions. However, because grains 
are not so strongly coupled with magnetic fields as ions due to their large masses, complete 
freezing is not a good approximation in most situations. More accurate treatment is required for 



their motion ( Elmegreen 1979 ). Moreover, the size and the mass of grains are distributed in wide 
ranges (e.g., Mathis, Rumpl, & Nordsieck 1977, referred to as MRN hereafter), and the degree of 
freezing depends sensitively on their masses. More accurate treatment is also necessary for ions 
because even ions are not well frozen at high densities. Nakano (1984) and Nakano & Umebayashi 
(1986, referred to as NU86| hereafter) formulated a method of describing the drift of magnetic 



fields in clouds containing any kinds of charged particles which are coupled with magnetic fields 
at arbitrary strengths. 

Using this formalism we investigated the time scale of magnetic flux loss, ts, from a major 
part of a cloud (core) as a function of the mean density of the cloud (core) ( [Nakano 1984 ; NU86| ; 



Umebayashi fc Nakano 1990| ; Nishi, Nakano, & Umebayashi 1991, referred to as NNU91 hereafter 



We have found that there is a critical value n& ec ~ 10 11 hydrogen nuclei per cm 3 for the density 
of the cloud, «h, at which tg is equal to the free-fall time tg of the cloud (core); ts < tg holds 
only at kh > "^deo and > at nn € ^dec- We have called this the decoupling density because 
extensive flux loss occurs only at nn ^ ^dec- 

Our previous results show that at 10 3 cm -3 <^ nn <^ 10 cm" 3 , is 10 to 10 2 times iff for 
B = B CT and a relation ts oc B~ 2 approximately holds at least for B cr > B ^ 0.1B cr . The relation 
oc B~ 2 is characteristic to ambipolar diffusion which occurs when the dominant charged 



particles are well frozen to magnetic fields (see § 2.2 for details). As the density increases at 
«h ^> 10 7 cm -3 , tB becomes more and more insensitive to B, and finally at nn ~ n& ec , becomes 
almost independent of B, reminiscent of Ohmic dissipation. Recently Desch & Mouschovias 
(2001) wrote that ambipolar diffusion was the dominant process even at densities several orders 
of magnitude higher than n^ cc . These results may cause confusion. Besides, the dependence of 
ts on nn and B described above is rather complicated. It would be necessary to clarify what is 
going on especially at high densities. Furthermore, although it was pointed out that grains play 
an important role in the process of magnetic flux loss, it does not seem to be widely recognized 
how important they are. 

The purpose of this paper is to clarify detailed mechanisms operating in the loss of magnetic 
flux in molecular clouds especially at high densities and to show how the grains behave. In § 2 we 
summarize some of the formulae obtained by Nakano (1984) and NU86| , which will be used in this 



paper. In § 3 we show numerical results and analyze in detail the processes going on especially at 
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high densities. Discussion is made in § 4, and summary is given in § 5. In Appendix [A| we show 



another method of the formulation than that of Nakano (1984) and NU86 



2. Drift of Charged Particles and Magnetic Fields 
2.1. Formulae 



We summarize some of the formulae obtained by Nakano (1984) and NU86| , which will be used 



in this paper. We consider a cloud which is composed mainly of neutral molecules and atoms but 
contains a slight amount of charged particles such as various atomic and molecular ions, electrons, 
and grains of various sizes and charges. 

Because each kind of charged particles are scarce, we can neglect in its equation of motion 
the pressure force, the gravity, and the inertia term compared with the Lorentz force and the 
frictional force exerted by the neutrals. With this approximation we obtain the drift velocity 
(velocity of the guiding center) v\ of an arbitrary particle A of mass m\ and electric charge q\ 
relative to the neutrals. We adopt the local Cartesian coordinate system with the z-axis parallel 
to the local magnetic field vector B and the x-axis along the magnetic force j X B/c, where 
j = (c/47r)V X B is the electric current density and c is the light velocity. The components of the 
drift velocity v\ are given by 

1 + (t a u; a ) 2 A{ + A\ \ t\lo\ J c 

Here, t\ is the viscous damping time of the motion of particle A in the sea of the neutrals, whose 
expressions are given, e.g., by Nakano (1984) for various particles, 

L0 X = (5) 

m\c 

is the gyrofrequency (defined to be negative for negatively charged particles), 
4 PvT v wl n u q u T u uj u ( B\ 2 

Al = Vi + (w = "4ri + (^) 2 = w CTP ' (6) 

^ p v u v B ^ n v q v ( B\ 2 

M = \ 1 + {r v u v f =7^1 + (r^f = UJ aH ' (?) 

where p u and n v = p v /m v are the mass and number densities, respectively, of particle v, and 
crp and <7h are Pedersen and Hall conductivities, respectively (see Appendix |A]). Summation in 
equations (||) and (0) is for all kinds of charged particles; particles having different values of t v or 
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u) v are different kinds. The component v\ z is not necessary because the drift along field lines has 
no effect on the magnetic flux loss. 

The drift velocity of magnetic fields, vb, is defined as the velocity relative to the neutrals of 
each point on an arbitrary closed contour in the cloud with which the magnetic flux through the 
contour is conserved. This requires 

E ± + -v B xB = 0, (8) 

c 

where E± is the component of the electric field E perpendicular to B in the frame moving with 
the neutrals. Using the relation between v\ and E, or the equation of motion, we obtain 

M 1 , • 

VB * = A\ + Al c |j X [ ' (9) 

M l , . „, 
VB y = ~ A\ + Ai ~c^ 3 x (10) 



The rate of magnetic flux loss is determined by vbx alone ( [NTJ86 ). This is self-evident for 



axisymmetric clouds because the local y-axis is in the azimuthal direction. Therefore, the subscript 
x of vbx an d v\ x may be omitted in the following. 

Using equations @ and (||) we can calculate the electric current density in the (x, y) plane, 
J = J2\ n \QxV\- With the electrical neutrality relation J2x Pa^a = 0, we can easily show that 
J x = 0, which is consistent with the definition of the x-axis, and that J y is equal to the component 
perpendicular to B of j which appears in equations (|3|) and (|j). Thus the formulation has been 
consistently done. Consistency of equations (||) and ( |l0| ) with equations (|3|) and (H) can be 
confirmed by considering the motion of a test particle which is completely frozen to magnetic 
fields and drifts with magnetic fields. Taking a limit Itxcj^I - > oo reduces equations (^) and to 
equations @ and (|i~0D, respectively, because A± and A2 are not affected by the test particle. 

In this formalism we have neglected the effect of charge fluctuation of grains caused by sticking 
of ions and electrons on their motion because Nakano &; Umebayashi (1980) and NU86| had a 



rough estimation that this effect is small. More elaborate discussion on this effect will be given 
in §^]^. We have also neglected the collision between charged particles, whose effect is negligibly 
small compared with the effect of their collision with the neutrals as shown in Appendix [B[ 

With the approximation made on the motion of charged particles, the magnetic force on 
charged particles must balance with the frictional force exerted on them by the neutrals, or by 
components 

E^ = -iixs|, (id 



yPxVxy =Q (12) 

A T * 

It is easy to confirm that v\ given by equations (||) and (Q) satisfies these equations. 
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One may simply think that the quantity |tawa| would characterize the degree of freezing of 
particle A to magnetic fields. However, we have from equations (H) and <M) 



V\x 
VBx 



1 + 



A 2 \ 



1 + (t a u;a) 2 V Mt\u\ 



(13) 



When |ta^a| S> 1, the degree of freezing, or the relative drift velocity, is certainly determined 
by t\uj\ alone as v\ x /vb x ~ (t\uj\) 2 /[I + (ta^a) 2 ] because IA2I/-A1 is at most several at least 
in the ranges of density and field strength covered by this paper. However, when |taWa| is n °t 
much larger than 1, the relative drift velocity may be greatly affected by the second term in the 
parentheses of equation (|i~3|), or by other charged particles. We shall show some examples of this 
effect in 3.4 and 3.5. 



Making use of equations (@), (0), and (10), NU86| showed that the electric field 
Eq = E — u n X B/c in the frame wherein the neutrals move with a velocity u n is given 
by 



E = --u n xB + —j + (3j X B-i{j X B) X B, 

c an 



where 



Cll 



D 



n u q p T u uj u 



is the electric conductivity parallel to magnetic field lines, and 



ft 



e-4 



B An 



c 2 A\ + Al> 
At 1 



(14) 



(15) 



(16) 



(17) 



c 2 A\ + A\ B 2 g\\ ' 

The work j • E gives the dissipation rate of magnetic energy per unit volume. The first term on 
the right-hand side of equation (14) causes amplification of magnetic fields by fluid motion. The 



second term leads to Ohmic dissipation. The third term yields no work. The last term gives rise to 
the dissipation in excess of Ohmic dissipation; ambipolar diffusion in excess of Ohmic dissipation 
is given by this term when the dominant charged particles are well frozen to magnetic fields. Thus 



the ratio of the last term to the second term in equation (14) 

criicrp 



D 



4 



4 



1 



(18) 



gives the excess dissipation relative to Ohmic dissipation flNU86p . 

In Appendix [A] we show another method of obtaining these formulae than that adopted by 
Nakano (1984) and |NU86. 
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2.2. Limiting Cases 

Because the drift velocity vb x is expressed in terms of the complicated quantities A\ and A2, 
it would be worthwhile to show how magnetic fields drift in some limiting cases. 

When the dominant charged particles are well frozen to magnetic fields, or \t v <jj v \ ^> 1, 
equations (||) and (0) give A\ « J2uPv/ T v 3> I-A2I ~ I J2u Pv/{ T v u v)\j and thus these particles drift 
at almost the same velocity irrespective of their charges as seen from equation (|3|). These relations 
also reduce equation @ to 

v Bx J2-*-\3XB\- (19) 

u T u C 

This can also be obtained by replacing v\ x with vbx in equation (|j"l|). Equation (|H]) means that 
the charged particles drift together with magnetic fields with the terminal velocity at which the 
magnetic force balances with the frictional force exerted by the neutrals. This is the generalization 
of ambipolar diffusion investigated first by Mestel & Spitzer (1956), who considered only ions as 
the transmitter of the magnetic force to the neutrals. Because the magnetic force balances with 
the gravitational force perpendicular to field lines when B = B CT , the mean magnetic force in an 
oblate cloud with B < B cr can be given by 

1 , / B \ 2 7T 

- I j X B| ft* (-^— ) 2 (») 

where p and U are the mean density and the mean column density along field lines, respectively, 
of the cloud. From equations ([II]) and (p0|) we have Vbx ^ B 2 , and then the magnetic flux loss 
time from the major part of the cloud, defined as the time required to drift a length scale L of 
magnetic fields, t B « L/v Bx oc B~ 2 . This is characteristic to ambipolar diffusion in magnetically 
supercritical clouds. 

When ions are the major contributors to the left-hand side of equation (p^), whose mass 
m\ is much larger than that of an H2 molecule and whose viscous damping time is given by 
T; ~ m\/(p{av)\) with the collision rate coefficient {av)\ ~ 1.5 x lCP 9 cm 3 s -1 (Nakano 1984), we 
have 



irGp V B 

for nearly spherical clouds with p ~ U/(2L) insofar as B < B cr . This is a well-known expression 
of the ambipolar diffusion time for magnetically supercritical clouds, which is proportional to the 
ionization fraction Wi/nn- 

If there exists a magnetically subcritical cloud core (B > B CT ) in equilibrium owing to external 
magnetic fields which suppress its expansion, the magnetic force therein almost balances with the 
gravity, or \j X B\/c ~ TrGSp/2 though this must be significantly smaller than B 2 /(AirR), where 
R is the core radius. Substituting this relation into equation di9|), we find that vbx and t B are 
independent of B and equal to those at B = B CT as shown by Nakano (1998). 
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Another limiting case is where the dominant charged particles are not at all frozen to magnetic 
fields, or \t u lu u \ <C 1. In this case we have 

f B\ 2 

A 1 &Y,P v t v uI= (-J a h (22) 
A 2 ^J2 ~ fal) = " E P^l (23) 

V V 

Electrical neutrality derives the last expression of equation (E3) . Again we find A\ 3> | A 2 \ ■ Taking 
| j X B\/c ~ B 2 /(4ttL), we obtain the drift time of magnetic fields 

L 47T(T|| „ , s 

This is the well-known Ohmic dissipation time, which does not depend on the field strength. 
Thus our formalism contains ambipolar diffusion and Ohmic dissipation as limiting cases. 



3. Numerical Results 
3.1. Cloud Model 



We adopt the cloud model almost the same as |NNU91| . Because magnetized clouds can 
contract along field lines rather easily, the force balance may approximately hold between gravity 
and gas pressure along field lines even when they are contracting dynamically perpendicular to 



field lines (e.g., Scott fc Black 1980| ; Nakamura, Hanawa, & Nakano 1995). In such clouds the 



mean density p and the mean column density U along field lines have a relation 

S ss — H ) ^ 0.040 . H „ — — gem" 2 , (25) 

\irGpmu) V10 5 cm- 3 10K p J B v ; 

where T is the mean temperature of the cloud, k is the Boltzmann constant, \i is the mean 
molecular weight of the gas, and mu is the mass of a hydrogen atom. Equation (|25|) gives the 
column density of an isothermal disk in equilibrium if p is half the density at the midplane. If 



magnetic fields are weak and the isothermal (spherical) cloud is on the verge of collapse ( Ebert 



1955 ; Bonnor 1956 ), the mean column density is 0.9 times that of equation (f2q). The half-thickness 



of the cloud is given by Z m E/(2p), or 

kT \ 1/2 nno „/10 5 cm- 3 T 2.37 \ 1/2 



Both S and Z are independent of the cloud mass. 
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3.2. Densities of Charged Particles 

As the reaction scheme of determining densities of various charged particles, we adopt the 
same model as Umebayashi & Nakano (1990) and NNU91| though we have revised some of the 



rate coefficients according to Le Teuff, Millar, & Markwick (2000); the greatest change is for 
dissociative recombination of H^. In dense clouds sufficiently opaque to interstellar ultraviolet 
radiation, ions and free electrons are formed mainly by ionization of H2 molecules and He atoms 
by cosmic rays. After some reactions in the gas phase and at grain surface they finally recombine 
each other. We determine the densities of various charged particles assuming steady state for all 
the reactions. 

As cosmic rays go deep into a cloud, their intensity decreases by interaction with matter, and 
the ionization rate decreases exponentially as £ = Co exp(— x/xq)j where Co is the ionization rate 
of an H2 molecule at the cloud surface, x is the depth in column density from the cloud surface, 
and xo ~ 96 gem -2 is the attenuation length of the ionization rate (Umebayashi & Nakano 1981). 
We take \ = because we are interested in the mean densities of various particles in the cloud. 
We also take into account the ionization by radioactive elements contained in the cloud at a rate 



6.9 x 10 23 s 1 ( Umebayashi fc Nakano 1981| ), which is important only at nn ;> 10 15 cm 



3 



As charged particles we consider electrons e~, atomic ions H + , He + , and C + , metallic ions 
such as Mg + , Si + , and Fe + , which we denote as M + collectively, H3" , molecular ions other than 
(typically HCO + ), which we denote as m + , and charged grains. As the mean masses of M + and 
m + we take 34mn and 29mn, respectively. We separate from m + because the large difference 
in their masses causes considerable differences in t\ and u)\. We adopt the MRN size distribution 
of grains given by 

diTt 

= Cn H a" 3 - 5 , a min < a < a max , (27) 

da 



where a is the grain radius and C 1.5 X 10 25 cm 2 - 5 ( praine fc Lee 198"4| ; Mathis 1986|) . We take 



a min = 5 nm and a max = 250 nm. As for the charge states of grains we consider neutral, ±e, ±2e, 
and ±3e, where e is the elementary electric charge. We divide the grain radius into 40 bins of 
equal logarithmic width Z\loga = 0.0425. This must be accurate enough for our purpose because 
we have found little difference in numerical results compared with the case of 20 bins adopted 
by |NNU91pl though we have found some difference from the case of 5 bins adopted by Desch & 
Mouschovias (2001). 

Figure 1 shows abundances of various particles as functions of the cloud density «h for the 
standard case Co = 1 x 10 _17 s _1 ( |Spitzer fc Tomasko 1968 ). Ions and electrons are dominant 



Desch & Mouschovias (2001) criticized NNU91 in their §2.1 claiming (i) ignorance of the size-dependence of 



grain charges, and (ii) ignorance of collisions between grains of different sizes which is important to neutralization 



of grain charges at high densities. In reality, however, NNU91 did not make these ignorances. As for (i) see their 
Appendix A, where the equations on grain charges have clear dependences on the grain radius, 
anything that means (ii) . 



NNU91 did not write 
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charged particles at wh <^ 10 6 cm 3 , ions and g (grains of charge — e) are dominant at nu between 
10 7 and 10 9 cm 3 , and g and g + are major at nu ^ 10 10 cm 3 . To keep electrical neutrality 
at very high densities nu ^ 10 10 cm -3 where ions and electrons are no longer dominant charged 
particles because of their efficient recombination, n.(g~) ~ n(g + ) must hold. Nakano (1984) 
showed that to assure this relation the number flux of thermal electrons must be (Si/Se) 1 ^ 2 times 
that of thermal ions, or 

n e ( Si m e \ l l 2 , . 

m \bcmj 

where S[ and S c are sticking probabilities of ions and electrons, respectively, when they collide 
neutral grains. Although Nakano (1984) ignored electric polarization of grains induced by 
approaching charged particles, which has an effect of enhancing the collision rate ( praine &z Sutin 



1987 ), this effect cancels out in the derivation of equation (|28|). As seen in Figure 1, metallic ions 
are dominant among various ions, or n\ ~ n(M + ), at all densities. In numerical calculation we 
took S e = 0.6 and 5; = 1 QUmebayashi fc Nakano 1980| ; |NNU91|) . This is why n e 6 x 10" 3 ra(M + ) 



holds at nn ^ 10 10 cm -3 . Equation (pg|) also yields n(g ) « n(g ++ ). Decrease of charged grains 
with nn at very high densities is due to neutralization of grains by mutual collisions. Cosmic rays 
are significantly attenuated at nu ^ 10 12 cm~ 3 . 

Figure 2 shows the charge state of grains as a function of the grain radius at several cloud 
densities. At low densities nu <; 10 6 cm -3 grains of charge — e are dominant and 1/4 to 1/30 of 
grains are neutral depending on their radius. At nu ~ 10 8 cm" 3 neutral grains are as abundant 
as — e grains and grains of other charges are much less at all radii. At very high densities 
nu ^ 10 10 cm -3 neutral grains are most abundant, and +e grains are as abundant as — e grains at 
all radii. 



3.3. Magnetic Flux Loss Time 

The mean magnetic force in an oblate cloud with B < B CT satisfies equation (^TJ) whether or 
not the force balance holds along field lines. The magnetic flux loss time ts can be given by the 
time for field lines to drift the length scale L of the cloud. From equations (^) and (|20| ) we have 

t B « « -g 1 2 29 

vbx \ B J Ai TrGUp 

First we consider a cloud wherein the force balance approximately holds between gravity and 
pressure along field lines and whose column density £ and half-thickness Z are given by equations 
( |25| ) and ((2?J), respectively. From equations @ and (^) we have 

/vr/cT/A 1 / 2 ,( T nu 2.37X 1 / 2 . . 

5 ^ 6.4 x 10~ 5 — - . a - G. (30) 

V iimu J V 10 K 10 5 cm~ 3 a J y ' 

The quasistatic contraction of clouds induced by the drift of magnetic fields is highly 
nonhomologous; only the densest central part of the cloud contracts leaving the outer part almost 



- 11 - 



unchanged ( [Nakano 197S , 1982; Lizano &z Shu 19~89^ ). The time scale of such contraction is the 



drift time of magnetic fields in the central part of the cloud whose length scale L is nearly equal 
to the thickness of the cloud, Z, or 

[B CT \ 2 A 2 + A 2 1 . , 

tB * (if) tt^' (31> 

Disk-like clouds in runaway collapse with field lines perpendicular to the disk layers have nearly 
uniform cores whose column density along field lines and size across them are given approximately 
by equations (|25|) and (pq) , respectively ( Nakamura et al. 1995| ). The flux loss time of these cores 



is also given by equation (|3l|). This equation also holds for nearly spherical clouds wherein the 
pressure force almost balances with the gravity (though B 2 must be significantly smaller than B 2 T ) 
because the cloud radius is nearly equal to the half-thickness along field lines, Z. 



As mentioned above, equation (|29|) for ts holds for dynamically contracting clouds even if 
the force balance does not hold along field lines. How about equation ([3l])? If the cloud is nearly 
spherical, equation ( |3l| ) can also be applied because with S/(2L) as p equation ( |29"|) is reduced to 
equation (|3l|). However, equation (||) should be used for B cr instead of equation (|30[), which holds 
only for clouds in force balance along field lines. 



The time scale ts given by equation ( |31| ) depends on the density and the magnetic field 
strength of the cloud, but does not depend on the cloud mass. The gas temperature T affects ts 
through densities of charged particles, t\, and B cr . We take T = 10 K in this paper. 

Figure 3 shows ie for the cases of B = B cr (solid lines) and B = 0.1B cr (dashed lines) for 
clouds in force balance along field lines, or with B cr given by equation (^). The dot-dashed lines 
show the Ohmic dissipation time t Q d, which is obtained by taking a limit B — > in equation (|3~ll), 



or by setting L ~ Z in equation (24). A relation tB = t Q d/(l + D) holds, where D is the excess 
dissipation relative to Ohmic dissipation given by equation (|l^). In addition to the standard case 
Co = 1 x 10 -17 s -1 (thick lines), we also show the case of Co = 1 x 10 _16 s _1 (thin lines) because 
there are some suggestions from observations of molecular ions that Co might be significantly 
larger than the standard value in some clouds (de Boisanger, Helmich, & van Dishoeck 1996; 



Caselli et al. 1998] ). At bh <^ 10 7 cm 3 , ts is 10 to 10 2 times % for the case of B = B CI , and ts 



for the case of B = 0.1-B cr is about 10 2 times larger than that for B = B CT . At njj <; 10 8 cm -3 , 
the ratio ts/tg decreases as the density increases, and at last ts = tg is attained at some density, 
which we denote n^ec and call the decoupling density as in our previous work. For the standard 
case Co = 1 x 10 _17 s _1 we find nd ec ~ 2 x 10 11 and 3 x 10 11 cm~ 3 for B = B CT and 0.1-B cr , 
respectively. As shown by |NNU91 , ts and n& ec do not depend sensitively on various parameters 



such as fractions of heavy elements remaining in the gas phase and on the details of the grain 
model (e.g., a m \ a ). As seen from Figure 3, ts and nd ec are not very sensitive to Co- 



Ciolek &: Mouschovias (1993) criticized our previous work (e.g., NU86 ; Umebayashi fc Nakano 



1990; NNU91) in their § 1.1 claiming that comparison of ts with ts and comparison of vb with 



the free-fall velocity ur were meaningless because magnetically subcritical clouds did not contract 
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freely. However, because the free-fall time is one of the fundamental time scales of clouds, we can 
find out by comparing ts with ts (or vb with ug) whether the magnetic flux is lost effectively in 
dynamically contracting clouds, and how slowly the clouds in quasi-equilibrium contract induced 
by the drift of magnetic fields without detailed simulations which have been done by many authors 
(e.g., [Nakano 1979| ; |Lizano k Shu 198$ piolek & Mouschovias T994[ ). According to Ciolek & 
Mouschovias (1993), we concluded erroneously that ambipolar diffusion (or drift of magnetic 
fields) was inefficient at tt-h <C nd ec because ts S> iff- Our conclusion in the previous and this 
papers is that extensive flux loss occurs only at tt-h <; ^dec- If there exist highly magnetically 
subcritical clouds with ^ 3> ^ C n they might lose magnetic fluxes extensively even at «h <C n^ec 
contracting quasistatically, though only down to ~ <£ cr , as shown by numerical simulations of 
Ciolek & Mouschovias (1994). However, clouds or cloud cores with <P 3> <& CT cannot exist flNakanc 



1998). Moreover, as mentioned in § 1, cloud cores must decrease their magnetic fluxes down to 



10 -3 # cr or below in the process of star formation. Cloud cores with somewhat smaller than 
'Per can begin dynamical contraction rather easily (Nakano 1998). Extensive flux loss down to 
^ -C ^cr does not occur at nn <C ndec because ts is much larger than the dynamical time scale. 



3.4. Dependence on Magnetic Field Strength 

We try to approximate the dependence of ts on the field strength B by a power law ts oc 5~ 7 . 
Figure 4 (top) shows the power index 7 obtained by comparing ts for the two cases B = B CT 
and 0.1i? cr for the standard case £o = 1 x 10 s _1 . At nn <^ 10 'cm -3 , deviation of 7 from 2 
is small, or ts oc B~ 2 approximately holds, at least for B CT > B ;> 0.1B cr , in good agreement 



with ambipolar diffusion described in § |2.2| . This is because ions and smallest grains, which are 



major contributors to transmitting the magnetic force to the neutrals, are relatively well frozen to 



magnetic fields in this density range. We shall discuss the small deviation of 7 from 2 in § 3.5. 



At ran > 10 7 cm -3 , 7 decreases steeply as the density increases, and finally settles down to 
0. The situation 7 <C 1 suggests that Ohmic dissipation contributes at least significantly to 
magnetic flux loss. However, although 7 = 0.13 at nu = 2 x 10 11 cm -3 , the flux loss occurs much 
faster than by Ohmic dissipation; ts ~ 0.1i o d as shown in Figure 3, or D 10. The steep decrease 
of 7 with density and this behavior of magnetic flux loss at nn ~ ^dec can be understood by 
checking the motion of some typical charged particles. 

Figure 4 (bottom) shows the drift velocity relative to that of magnetic fields, v\Jvb, given 
by equation (|l^), and t\uj\ of the dominant ions M + for the two cases B = B cr (solid line) and 
0.1-B cr (dashed line) with (0 = 1 x 10~ 17 s -1 . Here we omit the subscript x to the velocities. At 
nn = 1 x 10 9 cm -3 , for example, metallic ions drift almost with magnetic fields with v\/vb ~ 0.96 
for B = B cr because they are well frozen to magnetic fields with ~ 32. For B = 0.1B CT , 
however, we find v\)vb ~ 0.26 and t\uj\ ~ 3.2. This value of v\/vb is much smaller than a naive 
estimation v\/vb ~ (rjo;i) 2 /[l + (t\uj{) 2 ] ~ 0.91 for this value of T\u>i. This large deviation is caused 
by the second term in the parentheses of equation (|1|), or by interaction with other charged 
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particles, especially grains. Even for B = B CI deviation of V{ from vb is caused mostly by this 
term: a naive estimation gives v\/vb ~ (rio;i) 2 /[l + (rjCJi) 2 ] 0.999. As B decreases, v\/vb 
decreases markedly. 

The situation is the same for grains because |T g u;g| -C t x lo\. Figure 5 shows the drift velocity 
relative to that of magnetic fields, v g /vB, and |T g u; g | of grains of charge — e as functions of the 
grain radius at several densities for £o = 1 x 10 _17 s -1 . At n H = 10 9 cm 3 and B = B cv , grains of 
a = 10 nm with electric charge — e and +e, for example, have relative drift velocities v g /vB = 0.41 
and —0.22, respectively. These particles have |T g u; g | = 0.32, and then the naive estimation gives 
v g /vB ~ (T g u g ) 2 /[1 + (r g u>g) 2 ] 0.093 for both charges, much smaller than the actual values; the 
actual velocity of +e grains has even an opposite sign. As B decreases, \v g /vs\ decreases markedly 
for all radii as shown in Figure 5. 

At «h > 10 7 cm~ 3 ions are not very well frozen to magnetic fields even for B ~ B cr . Therefore, 
as B decreases, the degree of freezing and v-Jvb decrease markedly. The situation is the same for 
grains which are much less tightly coupled with magnetic fields. Decrease of v\/vb and \v g /vs\ has 
an effect of decreasing ts if V{ and v g are fixed. On the other hand, as B decreases, the magnetic 
force, or the driving force of the drift, decreases, which has an effect of decreasing the drift velocity 
v\ and \v g \ as seen from equation (|TT|), and then has an effect of increasing ts if v\/vb and v g /vs 
are fixed. These opposite effects make ts less sensitive to B than at ran < 10 cm" 3 , where ions 
and smallest grains, the dominant charged particles, are relatively well frozen and the decrease of 
B has little effect on their v\/vb- As the density increases, ts becomes more and more insensitive 
to B because |t^u;a| decreases with density for all charged particles even though ts is significantly 
smaller than t Q( j. This is why 7 decreases steeply with density at ran > 1 x 10 7 cm" 3 . 



3.5. Contribution of Grains 

Even at kh < 1 x 10 7 cm -3 , 7 shows some deviation from 2 as seen in Figure 4 (top). Because 
ions are strongly coupled with magnetic fields and deviation of Vi from vb is very small at these 
densities as shown in Figure 4 (bottom), this deviation suggests that transmission of the magnetic 
force to the neutrals is significantly contributed by grains, which are not very strongly coupled 
with magnetic fields even at these low densities as shown in Figure 5. 

To confirm this we calculate the frictional force exerted by each kind of charged particles on 
the neutrals, which is given by each term of equation (|TT|) . We shall not consider the component 
of the frictional force perpendicular to both B and j X B because summation of this component 



for all charged particles vanishes as shown by equation (12). Figure 6 shows the frictional forces of 
ions and grains relative to the total frictional force, which is equal to | j X B\/c per unit volume, 
as functions of the cloud density for the two cases of the field strength with £0 = 1 x 10" 17 s~ 1 . 
Ions contribute more than grains only at ran <^ 10 4 cm~ 3 , where the ionization fraction is relatively 
high. At n H ^ 10 6 cm 3 grains contribute more than about 90% of the frictional force at 
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least for B cr > B ;> 0.1i? cr . The frictional force exerted by electrons is more than 3 orders of 
magnitude smaller than that by ions mainly because of their small drift momentum. For the case 
of Co = 1 x 10" 16 s" 1 these densities are somewhat higher; e.g., ions contribute more than grains 
at wh ^ 10 5 cm -3 because of higher ion densities. 

Figures 7 and 8 show the frictional force exerted by grains on the neutrals as a function 
of their radius a and charge q g at several cloud densities for the two cases of the field strength 
B = B CI and 0.1-B cr with £o = 1 x 10 s . As well as ions, negatively charged grains exert 
frictional force, or drift, in the direction of the magnetic force irrespective of their radius at 
n H ^ 10 14 (-B/l? cr ) 2 cm -3 at least for B ^ 0.1B cr . Positively charged grains behave differently. 
At low densities they contribute in the same direction as negatively charged grains irrespective 
of their radius. However, at some density the largest grains begin to contribute, or drift, in the 
opposite direction, and the size range of such grains expands to smaller radii as hh increases; 
frictional forces in these ranges are shown by dotted lines in Figures 7 and 8. For example, grains 
of charge +e with a > 160 nm drift in the opposite direction at rin = 10 5 cm" 3 for B = B cr , and 
those with a > 93 nm do at ran = 10 4 cm" 3 for B = 0.1-B cr though the dotted line is outside this 
panel. Finally at nn ~ 2 x 10 9 and 9 x 10 6 cm" 3 for B = B CT and 0.1-B cr , respectively, even the 
smallest grains (a ~ 5nm) with charge +e begin to drift opposite to the magnetic force. These 
are because Ai < (or cjh < 0) in most of the density region covered by this paper (see Figure 
9; A\ > by definition), and thus positively charged particles with small t\lu\ have negative drift 
velocities as seen from equation (13). Because T g u g oc a" 2 n H 1 ^ 2 (B / B CT ) , the radius range in which 



positively charged grains have negative drift velocities expands to smaller a as ran increases for a 



fixed B/B cr as long as A2 < 0. We shall discuss the physical reason of this phenomenon in § |4.1| . 

At high densities, not only ts deviates greatly from the ts oc B~ 2 relation, but also grains, 
the main contributors to the frictional force, of opposite charges drift in opposite directions. This 
contradicts the literal meaning of ambipolar diffusion, a term originally used in plasma physics for 
a quite different phenomenon in which particles of opposite charges diffuse as one unit (e.g., Cap| 
1976 ; see also Appendix |C|). 



The great contribution of grains to the frictional force comes from their large contribution 
to Ai and A2, or op and cjh, though <tm is mainly contributed by electrons in most of the 
density region in this paper. This difference in the contribution is due to large differences 
in |7au;a| an d in the abundance n\/nn among charged particles. For metallic ions we 
have TitJi xi 320(n H /10 7 cuT 3 )- 1 / 2 (B / B^) . At T ps 10K we find |r e cj c |/riCJi ps 4.8 x 10 3 and 
\T g U!g\/TiUJi 1.0 x 10" 2 (a/10nm)" 2 for grains of charge q g = dbe. For the case of Co = 1 x 10" 17 s , 
a\\ is mainly contributed by electrons at ran ^ 10 13 cm -3 and by grains at higher densities as shown 
in Figure 9. Contribution of ions to cry is minor at all densities. To crp or Ai, ions contribute 
more than half only at bh ^ 10 cm -3 , and grains contribute more than 80% at re H £ 10 5 cm 3 for 
B = B CT (Figure 9), though for B = 0.1B cr ions also contribute 50 — 80% at bh between 3 x 10 8 and 
6 x 10 10 cm -3 . At rtn 10 11 cm" 3 , grains overwhelm the other particles. Contribution of electrons 
to cp is less than 10 -3 because of large |r e u;e| at low densities and because of low abundance at 
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high densities [see equation (||) or flAq)]. Situation is complicated for cjh or Because the terms 
of some given particles have quite different forms even with opposite signs between the equivalent 
equations (^) and flA6|), we cannot uniquely tell how much each kind of particles contributes to <th- 
Moreover, particles of opposite charges contribute oppositely and significant cancellation occurs in 
some situations for either expression of cjh- We would make a serious mistake if we neglect the 
terms of grains in <th except at high densities where both TiLO\ <J 10 and ri\ 3> n e hold. 

The sign of cth changes at some density as shown in Figure 9. With an expression of <th 
accurate at very high densities, we find that this occurs when t\uj\ ~ (ric/re;) 1 / 2 . As long as 
B ^ 0. 01i? cr , this gives T\U!\ ~ 0.08, which corresponds to the density bh ~ 1.6 x 10 14 (i?/-B cr ) 2 cm~ 3 . 
Because of the change of the sign, grains change the directions of drift at densities somewhat 



higher than this as seen from equation (13). However, the total frictional force they exert on the 
neutrals is always in the direction of magnetic force because the mean drift velocity of — e and +e 
grains is in the direction of magnetic force as seen from equation (J3T) or (|l3|). 



Desch & Mouschovias (2001) ignored the contribution of grains to a-p (<7j_ by their notation) 
and <7h as seen from their equations (25) and (26). Therefore, their equation (28), which was 
obtained by using these equations, does not correspond to the critical state at which Ohmic 
dissipation becomes important, or D ~ 1, contrary to their statement. Their results will be 



discussed as compared with ours in §|4,5 



Because the frictional force is contributed mainly by grains, is determined mainly by grains 
except at very low densities. Besides, the reactions at grain surface are important in determining 



the densities of various charged particles (e.g., Nakano 1984; NNU91). Thus, grains play a decisive 



role in the process of magnetic flux loss in molecular clouds. 



4. Discussion 
4.1. Drift of Grains 

Grains of opposite electric charges drift in opposite directions at high densities as shown in 



3.4 and 3.5. This is formally a result of the second term in the parentheses of equation (13). 



Here we shall give a more intuitive explanation of this phenomenon. 

In the existence of an electric field E a particle of electric charge q\ ^ gyrating around 
magnetic field linesf] drifts with a velocity v\ = cE X B/B 2 independent of its charge. If there is 
a force field /, which is independent of the particle charge, instead of the electric field, the particle 



3 One may think that gyration completely dissipates because the viscous damping times t\ are smaller than ts 
and iff by orders of magnitude. However, because of collisions with neutral molecules and atoms charged particles 
always have thermal motions, which would cause gyration, though a term "gyrating" may not be appropriate when 
\t\u\\ <; 1. 
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drifts with a velocity v\ = (c/q\) f X B / B 2 dependent on its charge. In the frame moving with 
the neutrals in a molecular cloud, there is an electric field satisfying equation (|8|), which causes a 
charge-independent drift motion (call this drift 1). Neutral molecules exert a frictional force on 
this drift motion, which is anti-parallel to the drift velocity independent of the electric charge. 
Therefore, this force causes another drift motion (drift 2) whose direction depends on the particle 
charge. Thus the drift velocity has both components dependent on and independent of the particle 
charge. Because drift 2 has opposite directions for particles of opposite charges, the frictional force 
on this motion yields the third drift motion whose direction is charge- independent, etc. 



All these can be taken into account in the equation for the steady motion ( Al) in Appendix [A], 



The direct solution of this equation is given by equations ( |A2| ) and (A3), whose components 
perpendicular to B are found to agree with equations (||) and @. When |t^cj^| <; 1, the frictional 
force can be comparable to or stronger than the electric force, and therefore the charge-dependent 
part of the drift velocity can be comparable to or greater than the charge-independent one. 

Because grains of opposite charges drift in opposite directions at high densities where grains 
are the major charged particles, one may anticipate that charge separation occurs in the j X B -, 
or x-direction. However, this is not the case. All charged particles should be taken into account. 



As confirmed in § 2.1, equation (^) guarantees that the x-component of the electric current density 
vanishes. 

Although the frictional force is contributed mostly by grains at high densities, their drift 
velocities are very small. For example, at the decoupling density ran ~ n^ec ~ 2 x 10 11 cm -3 for 
B = B CI with (o = 1 x 10 _17 s _1 even the smallest grains of a ~ 5nm with electric charge — e 
have a drift velocity as small as v g /vB ~ 0.02 as seen from Figure 5. Because extensive magnetic 
flux loss occurs only at ran ^ ridec and decrease of B / B CT by the flux loss makes v g /vs even 
smaller (Figure 5), grains are hardly lost from the cloud core even if the magnetic flux decreases 
by orders of magnitude as far as the cloud core is magnetically supercritical. Although Ciolek 
& Mouschovias (1994) write that the grain abundance decreases almost in proportion to the 
reduction factor of the magnetic flux in highly magnetically subcritical cloud cores, cloud cores 
can hardly be magnetically subcritical as discussed by Nakano (1998). Although the drift velocity 
of ions is not very small compared with that of magnetic fields even at «h ~ n^ cc for B ~ B cr 
(Figure 4), ions are quite minor constituents among the heavy elements even in the gas phase 
(Figure 1). Thus magnetic flux loss in cloud cores has little effect on the abundance of heavy 
elements in stars born therein. 



4.2. Charge Fluctuation of Grains 

Electric charges of grains change stochastically because ions and electrons sometimes stick to 
them. This change has some effect on their motion in electromagnetic fields. Nakano (1984) and 
NU86| neglected this effect because they had a rough estimation that it was small (see also [Nakano 
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& Umebayashi 1980). 

Kamaya & Nishi (2000) investigated this effect more elaborately for grains of charge — e. The 
probability of having some charge is proportional to the time span during which the grain has this 
charge. Therefore, as seen from Figure 2, the grains of charge — e change their charge mostly by 
going to the neutral state at most radii a and at most densities njj- The charge change between 
the — e and neutral states effectively decreases the viscous damping time r g of grains of charge — e 
by some factor C g > 1. This effect can be taken into account by replacing r g with r g /C g in our 
formalism summarized in §||. Kamaya & Nishi found that if is the dominant ion, C g for grains 
of a = 10 nm is 1.3 at bh ~ 10 3 cm -3 and approaches 1 as nn increases, and that larger grains 
have C g closer to 1. 

Because metallic ions and molecular ions other than are the dominant ions as shown 
in Figure 1, deviations of C g from 1 are about 3 times smaller than those obtained for by 
Kamaya & Nishi (2000) at all densities and at all grain radii; the deviation C g — 1 is proportional 

— 1/2 

to m i , where m\ is the mass of the dominant ions. Therefore, charge fluctuation of grains can 
hardly affect ts because C g can be as large as 1.1 only for smallest grains only at low densities 
nn ~ 10 3 cm -3 , where grains are minor contributors to the frictional force. Thus, our formalism 
can be applied to any situation with little error though Ciolek & Mouschovias (1993) pointed out 
our ignorance of this effect. 



4.3. Comparison with Other Work 

We compare our results with several previous works on ts, n^ cc , and the dissipation in excess 
of Ohmic dissipation given by D in equation (|l£ 



The decoupling densities n<j e c obtained in § |3.3| are about 4 times larger than that of |NNU91 

for the same grain model. The main cause for this difference is in B cr ; our B CT is y/2> times larger 

than theirs. At relatively low densities where major charged particles are well frozen to magnetic 

fields, difference in the coefficient of B CT has little effect on ts for a given B/B CT because these 

particles drift almost with magnetic fields and the drift velocity vbx is determined by the balance 

of the frictional force on them with (B/B cr ) 2 times the gravity as seen from equations ( |T^ ) and 

(|20|). At high densities, ts, and then n& ec , depend on the coefficient of B CT because the delay 

in the drift of major charged particles, which are not well frozen, depends on the field strength 

as discussed in § |3.4| . However, the general feature of ts (rough dependences on nn and B/B CT , 

existence of nd ec , etc.) is hardly affected by this coefficient. 
I I , 

NU86 investigated magnetic flux loss in spherical clouds, and found that magnetic decoupling 
occurs at nn « ra^ec ~ 5 x 10 11 cm~ 3 for a cloud of 1M and the decoupling density does not 
depend sensitively on the cloud mass. They also found that dissipation in excess of Ohmic one is 
not large, or D ~ 1, at nn ~ n dec for B « B cr . 
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In §|3.4| we have found D ~ 10 at wh ~ ndeo which is an order of magnitude larger than 



D obtained by NU86. We can list up some causes for this discrepancy; differences in the cloud 



model and in the grain model. NU86 assumed that clouds were spherical without force balance 
along field lines. They also assumed for simplicity that all grains had the same radius a = 100 nm. 
Another difference is in the electric conductivity an . For the collision of electrons with the neutrals, 



NU86| adopted the classical Langevin's cross sections in which electric polarization of the neutrals 



induced by an approaching electron has a great effect. These were widely adopted formerly 



(e.g., Nakano 1984 ; Mouschovias 1991 ) assuming the similarity with the collision of ions with the 
neutrals. However, we found that laboratory experiments show much smaller cross sections at low 
energies (payashi 198l| ), which are not much different from the geometrical cross sections. Sano et 



al. (2000) give the empirical formulae of the momentum transfer rate coefficients for the collision 
of electrons with H2 molecules and with He atoms, which were obtained by Umebayashi (1993, 
private communication) by fitting to the laboratory data. With these collision rate coefficients, 
which are about 80 times smaller than the previous ones at T ~ 10 K, the conductivity an, i di 
and D are much larger than the previous ones except at bh ^ 10 13 cnr 3 where an is mainly 
contributed by grains. On the other hand, ap is hardly affected by this change of the cross sections 
because contribution of electrons is very small even with the enhanced a e (Figure 9), and <th is 
also hardly affected except at very high densities where |<th| <C ap. Therefore, ts = i d/(l + D) is 
hardly affected by this change because i d oc an and 1 + D = a\\ap/(ap + a#). 

Desch & Mouschovias (2001) write that Ohmic dissipation becomes important, or D ~ 1 
is realized, when |r e u; e | ~ 6. However, |r e u;e| ~ 6 was obtained from their equation (28), 
TiCJi|r c u; c | ps n c /n;,0 which does not correspond to the state of D ps 1 because they ignored the 
contribution of grains to ap and cth (see § |3.5[ ). Their criterion |r e o; e | ps 6 for D ps 1 gives a density 
nn ps 6 x 10 17 (B j B CT ) 2 cm -3 irrespective of Co when B CT is given by equation fl30|), or the pressure 
force balances with the gravity along field lines. When contraction along field lines is insufficient 
and then the pressure force along field lines is weaker than the gravity, U is larger than given by 
equation ( |25|) and B cr is larger than equation (|3(i| ) for a given density. This means that \r e u) e \ is 
larger for given «h and B/B cr , and therefore |T e u;e| ~ 6 gives even higher densities than in the 
above case. 

Taking into account the contribution of grains to ap and cjh, we have found that D ~ 1, or 
*B ~ *od/2, is realized at «h ~ 1 x 10 13 cm 3 for Cq — 1 x 10 ^ s ^ almost independent of B j B cr 



4 Emphasizing the contrast with this equation, Desch & Mouschovias (2001) write that NU86 used a relation 
TiU>i\T c uj c \ ~ 1 for D ~ 1 assuming n c = rij even at densities as high as nu > "-dec- This is another misrepresentation 



of the fact in addition to their comments on 



NNU91 



cited in footnote 2. Although NU86 gave the expressions of D, 
by their equations (49) and (50), when ions and electrons are major contributors to A\ and A2 (or ap and an) (of 
course n c ~ m), and when grains of charge — e and +e are major contributors, respectively, as limiting cases, they 
used the general expression of D given by their equation (48) or ours ([L8|) in numerical calculation. For instance, to 
calculate D shown in their figure 2, they used the densities of charged particles shown in their figure 3, which clearly 
shows the deviation of n e from n, at high densities in agreement with our equation (|28|). 
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at least at 0.1 ^ B/B cr < 1. This can be confirmed by substituting the values of cry, ap, and <7h 
shown in Figure 9 into equation (^) and by noting that is almost independent of B/B CT at 
high densities as shown in Figure 3. This density is lower than Desch & Mouschovias' values by 
orders of magnitude. 

In reality, however, Ohmic dissipation becomes important at significantly lower densities. At 
n H ^, "dec the clouds contract dynamically adjusting B/B cr so as to keep ts ~ iff because too 
much decrease of B/B cr makes ts larger than iff as long as i d > tg. However, ts = t ^/(l + D) 
cannot be larger than t Q( j. Therefore, after i Q d ~ iff is attained at nn ~ 1 x 10 12 cm~ 3 (Figure 3), 
field dissipation proceeds by Ohmic dissipation faster than the contraction. 

Kamaya & Nishi (2000) investigated the motion of ions and grains in the process of magnetic 
flux loss using a simplified model. Assuming that all grains have the same radius and have an 
electric charge — e or 0, ions are well frozen to magnetic fields (riO>i 3> 1), and electrons are 
completely frozen (|r e o; e | — ► oo), they obtained the drift velocities of ions and grains of charge 
— e. If we apply equation (0) to their model, we easily obtain the drift velocities given by their 
equations (34) and (35) though r g should be replaced by T g /C g . Thus, Kamaya &; Nishi (2000) 
obtained some relatively transparent results due to their simplified model, which are consistent 
with our formalism, though they got some results different from Ciolek & Mouschovias (1993). 



5. Summary 

We analyzed the detailed processes operating in the drift of magnetic fields in molecular 
clouds taking into account charged grains with the MRN size distribution in addition to ions and 
electrons and the effect of partial freezing of these particles to magnetic fields using the formalism 



obtained by Nakano (1984) and NU86 



We found that to the frictional force, whereby the magnetic force is transmitted to neutral 
molecules, ions contribute more than half only at cloud densities rin <J 10 4 cm~ 3 , and grains 
contribute more than about 90% at n H £ 10 6 cm . Besides, the reactions at grain surface are 
important in determining the densities of various charged particles. Thus grains play a decisive 
role in the process of magnetic flux loss in molecular clouds. 



We confirmed the previous results (NNU91) on the magnetic flux loss time ts and the 



decoupling density ndeo at which i# is equal to the free-fall time iff and therefore magnetic fields 
are effectively decoupled from the gas; e.g., *b > f(f at wh < "-dec, and nd ec is almost independent 
of the magnetic field strength in the cloud and takes a value a few xlO 11 cm' 3 . We found that i# 
and nd ec are not very sensitive to the ionization rate by cosmic rays, £o- 

We investigated the dependence of i# on the field strength B. Approximating the relation 
by a power law ig oc B~' y , we found 7 ~ 2, characteristic to ambipolar diffusion, only at 
n u ^ 10 7 cm -3 where ions and smallest grains are pretty well frozen to magnetic fields. At 
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n H > 10 7 cm 3 , 7 decreases steeply as rin increases, and finally at nn ~ ^dec, 7 "C 1 is attained, 
reminiscent of Ohmic dissipation, though the flux loss occurs about 10 times faster than by Ohmic 
dissipation at njj ~ a few xlO 11 cnr 3 . Because even ions are not very well frozen to magnetic 
fields at nn > 10 7 cm~ 3 , ions and charged grains drift slower than the magnetic fields. Decrease of 
B has an effect of increasing this lag. This insufficient freezing makes ts less sensitive to B than 
a t «h ^ 10 7 cm -3 where major charged particles are well frozen and the decrease of B has little 
effect on their lag. This tendency is enhanced as the density increases, and at last ts becomes 
almost independent of B at tt-h ~ n^ ec . Ohmic dissipation is dominant only at nn ;> 10 12 cm -3 . 

We found that while ions and electrons drift in the direction of magnetic force at all densities, 
grains of opposite charges drift in opposite directions at high densities, where grains are major 
contributors to the frictional force, apart from the component of the drift velocities perpendicular 
to the magnetic force which yields no net frictional force. Although magnetic flux loss occurs 
significantly faster than by Ohmic dissipation even at nu ~ f^deo the operating process at high 
densities is quite different from ambipolar diffusion in which particles of opposite charges are 
supposed to drift as one unit. 

This work was supported in part by the Grant-in-Aid for Scientific Research of Priority 
Areas (A) (No. 10147101, 10147105) of the Ministry of Education, Culture, Sports, Science, and 
Technology of Japan. 



Another Method of Formulation 



Because the procedures of obtaining the formulae summarized in § 2.1 are rather complicated, 
it would be worthwhile to show another method of obtaining them. 



With the same approximation as adopted by Nakano (1984) and NU86 and described in § 2.1 
the mean motion of charged particle A obeys 

Je + 1 - Vx xb)-^ = o, (Al) 



in the frame moving with the neutrals (see § 24. for the notation). The solution of this equation is 
given by 

v n = q -^E h (A2) 

ii mx ii 

B 1 + (txujxY \ t\u\ B J 

where subscripts || and _L represent the components parallel and perpendicular to B, respectively. 

The electric current density j is obtained by summing up n\q\V\ for all kinds of charged 



particles using equations (|A2|) and (A3). Moving to a more general frame wherein the neutrals 
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move with a velocity u n and the electric field is given by Eq 
electrical neutrality relation, we obtain 



E — u n X B/c, and using the 



1 



B 



<T\\E \\ + cr P ( E 0± + -u n x B) + cjh— X ( £ _L + - u a X S ) , 



B 



1 



(A4) 



where try is the electric conductivity along magnetic field lines given by equation (^), and op and 
o"h are Pedersen and Hall conductivities, respectively, given by 



C \ 



BJ 



A 2 , 



(A5) 



(A6) 



where A\, A 2 , and ov are given by equations (||), ([?]), and (|15D, respectively. The last equality of 



equation (|A6|) can be found by using the electrical neutrality relation. Equation (|A4|) with (|15|), 
(|A5D , and (|A|) is a generalization of, e.g., Parks (1991) who considered only electrons and a single 
kind of ions as charged particles. 

The drift velocity of magnetic fields, vb, satisfies equation @. Making a vector product of 
equation ( A4 ) with B and eliminating Eq + u n X B/c = E by using equation (JsT) , we obtain 



1 



c 

This equation is solved for «b as 



<7p 



-j X B + 



0"H 



1 / • „x S 

2- 7 (j X S) X — 



H 



(A7) 



(A8) 



By taking the components of this equation we can confirm that equation (A8) is identical with 
equations (||) and (|To|) . 



(A9) 



Eliminating E± in equation ( |A3[) by using equation (||), we obtain 

(ta^a) 2 / 1 



1 + (ta^a) 2 

Elimination of vb by using equation ( |A8[ ) gives 



f A_L 



(ta^a) 5 



1 + (t a o;a) 2 \b) ap+a H 



l 



<Tp + 



TA^A 



-j xB 

c 



ap 
TA^A 



^IgX-O X B) 



(A10) 



It is easy to confirm that equation ( AlCf ) is identical with equations (||) and 
With some manipulation equation ( |A4j ) can be solved for Eq as 



1 1 . a H 
Eq = U n X B -\ J -\ — — 

C (Til Op + <7 H 



5 

2 + 



B 



1 \ -B 



al + a'i 



a\\ 



B 



x hx 



5 



(All) 



By comparing the coefficients we can confirm that equations ( |i~4| ) and (All) are identical. Desch 
& Mouschovias (2001) used an equation equivalent to ours (jil|) or ( [All ) though neglecting the 
dominant terms of <7p and <th as pointed out in 
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B. Effect of Collision between Charged Particles 

In §|2] and Appendix |A| we have neglected frictional forces on charged particles except 
those exerted by the neutrals. Except for electrons in some limited situations, this is a good 
approximation. We discuss the effect of collision between electrons and other charged particles on 



the collision time r e of electrons, which appears in the electric conductivity a e in equation (15). 

The 90° deflection time of an electron of velocity v e by collision with ions of electric charge 
q\ = e is given by 

Te -i o n — i' y m > 

OTrrii In A 

where ln/1 is a quantity determined by the cut-off of the impact parameter for the collision 
flSpitzer 1962] ). For v e we take the thermal velocity of electrons. The collision time of an electron 



with the neutrals is given by 



[ n(H 2 ) <r(e - H 2 ) + n(He) <r(e - He) ] v e ' 



(B2) 



where a(e-H2) ~ 6.6 x 10 _16 cm 2 and cr(e-He) ~ 2.9 x 10 _16 cm 2 are the collision cross sections 
of an electron with an H 2 molecule and an He atom, respectively, at T ~ 10 K ( [Sano et al. 2000| ). 
The collision time r e of an electron with ions or neutrals is given by 

111 

- = — + —. (B3) 



(B4) 



From equations ( [BID an d ( |B2| ) we have 



1 x 10- 10 J V 40 



for T ~ 10 K. Thus, the electric conductivity an is significantly affected by collision of electrons 
with ions only when n;/nn ^ 1 x 10 -10 , or at nn ^ 1 x 10 8 cm -3 as seen from Figure 1. 
Although charged grains also scatter electrons, their contribution does not affect this conclusion 
much because n(g + ) <^ n(g~) <^ n\ at bh ^ 1 x 10 8 cm -3 , [n(g~] + n(g + )]/nn ^ 1 x 10 -10 
at ^ 1 x 10 8 cm~ 3 , and lnyl for electron-grain collision is not much different from that for 
electron-ion collision. 

Similarly ions are scattered by charged grains. This effect is important compared with their 
collision with the neutrals only when [n(g~) +n(g + )]/nH ^ 2 x 10~ 9 , which is not realized as seen 
from Figure 1. Collision of ions with electrons has negligible effect on the motion of ions unless 
n e/ n H ^4x 10~ 7 . Collision of charged grains with ions is also much less effective than their 
collision with the neutrals. 

Although cry is greatly affected by the scattering of electrons by ions at kh <^ 10 8 cm~ 3 , up 
and cjh are hardly affected as can be confirmed in the following way. Even if the effect of collision 
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with ions is taken into account on r e at kh ^ 1 X 10 8 cm -3 , we find |r e u; e | 3> t\uj\ > 1 as far as 
nj/nn <C 10 -6 , which is satisfied at least in the density range covered by this paper. Therefore, 
contribution of electrons to A\ and A2 (or up and <th) is still much smaller than that of ions as can 
be confirmed from equations (||) and (0). Our results are therefore not affected by the correction 
on r c . On the other hand, we have used cni with the corrected r e in the estimation of the Ohmic 
dissipation time t Q( i « 10 15 yr at «h ~ 10 5 cm -3 in § 1; without this correction, t oc j at this density 
will be 10 2 times larger. 



C. On the Term "Ambipolar Diffusion" 

Ambipolar diffusion is a term originally used in plasma physics. At length scales larger than 
the Debye shielding length in plasmas, electrical neutrality is well established. Therefore, if there 
is a gradient in the electron density, they are transported with a flux proportional to their density 
gradient, and ions follow them because of electric forces, and vice versa. Thus, particles of both 
electric charges diffuse as one unit. This process is called ambipolar diffusion (e.g., |Cap 1976| ), 
whose literal meaning fits the process. 

In the drift of magnetic fields in molecular clouds investigated first by Mestel & Spitzer 
(1956), ions and electrons, which are well frozen to magnetic fields, drift with the same velocity 
as if they were one unit. Because of this seeming similarity this process got the name "ambipolar 
diffusion". However, the drift is driven by the magnetic force, not by the density gradient. Thus, 
this process is quite different from what the physical term "diffusion" means. 

In the original ambipolar diffusion the diffusion coefficient perpendicular to magnetic field 
lines is inversely proportional to both B 2 and the collision time of electrons with neutrals, r e , when 
t\uj\ S> 1 ( |Cap 19761) , an d thus the diffusion time is proportional to r e B 2 . On the other hand, the 
drift time of magnetic fields given by equation (|2l| ) in the same situation t\l>j\ 3> 1, ig oc t^ 1 B~ 2 , 
has the opposite dependences. This difference also shows that these processes are quite different 
from each other. 

Because of these differences the latter process was sometimes called "plasma drift" instead of 
ambipolar diffusion (e.g., [Spitzer 197S ; Nakano 1984 ; Mestel 1999 ). 



In this paper we have found that grains, the main contributors to the frictional force, of 
opposite charges drift in opposite directions except at relatively low densities. This is even 
contrary to the literal meaning of ambipolar diffusion. Some appropriate term is desired; plasma 
drift is much better than ambipolar diffusion. 
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Fig. 1. — Abundances of various particles, n(X)/nn, as functions of the density nu of the cloud by 
number of hydrogen nuclei. The solid lines are for ions, and the dotted line is for electrons. The 
dashed lines labeled g^ represent number densities relative to wh of grains of charge xe summed 
up over the radius. We have taken the ionization rate of an H2 molecule by cosmic rays outside 
the cloud, Co = 1 x 10~ 17 s _1 (standard case). We have assumed that 20% of C and O and 2% of 
metallic elements remain in the gas phase and the rest in grains. 
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Fig. 2. — The charge state distribution of grains as a function of the grain radius a at several cloud 
densities for the same parameters as in Fig. 1. Each panel is labeled with »h in cm -3 , and each 
line is labeled with the grain charge q g . 
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Fig. 3. — Time scales of magnetic flux loss for the clouds in which force balance approximately holds 
along field lines, or B CI is given by equation (|30|). The flux loss time is is shown for the two cases of 
field strength B = B CT (solid lines) and B = 0.1B CT (dashed lines). The Ohmic dissipation time t d 
is shown by the dot-dashed lines. Two cases of the ionization rate by cosmic rays, Co = 1 x 10 _ s _1 
(thick lines: standard case) and 1 x 10~ 16 s _1 (thin lines), are shown. The other parameters are 
the same as in Fig. 1. For comparison the free-fall time % = (3ir /2>2Gp) 1 / 2 is shown by the dotted 
line. 
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Fig. 4. — Dependence of the magnetic flux loss time ts on the field strength B and the drift velocity 
of ions for the standard case Co = 1 x 10~ 17 s~ 1 . We approximate ie by a power law tg oc B~ J . Top: 
The power index 7 obtained by comparing ts for the two cases B/B CT = 1 and 0.1 shown in Fig. 
3, 7 = — A log ts/A log B, as a function of the cloud density. For weaker magnetic fields we would 
obtain smaller values of 7. Bottom: The drift velocity in the direction of magnetic force relative 
to that of magnetic fields, Vi/vb, given by equation (|l~3|), and tiuji of metallic ions M + , dominant 
among various ions, for the two cases of field strength B = B cr (solid line) and B = 0.11? cr (dashed 
line). We have omitted the subscript x to the velocities. These quantities take almost the same 
values for molecular ions m + other than Hjj~, abundant next to M + , because their mean mass is 
not much different from that of M + . 
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Fig. 5. — The drift velocity relative to that of magnetic fields, v g /vB (top), given by equation (13), 
and |TgU;g| (bottom) for grains of charge — e as functions of the grain radius a at several densities 
for the same case as in Fig. 4. The two cases of field strength B = B cr (solid lines) and B = 0.1i? cr 
(dashed lines) are shown though the dashed lines for T g w g overlap with the solid lines except for 
lognn = 10 and 11 in cm~ 3 . The values of logriH are attached to the right ends of the lines for 
B = B CT and mostly to the left ends for B = 0.1-B cr . 
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Fig. 6. — Frictional forces exerted by ions and by grains on the neutrals relative to the total 
frictional force, which is equal to | j X B\/c per unit volume, as functions of the cloud density nn 
for the same case as in Fig. 4. Each kind of particles exerts the frictional force given by each 
term of equation (11). The frictional forces were summed up for all kinds of ions and for grains 
of all radii and all charges. Shown are the two cases of the field strength B = B cr {solid lines) 
and B = 0.1-B cr {dashed lines). The frictional force exerted by electrons is more than 3 orders of 
magnitude smaller than that by ions. 
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Fig. 7. — The frictional force exerted by grains on the neutrals as a function of their radius a and 
electric charge q g at several cloud densities for the same case as in Fig. 4 with B = B CI . Each 
curve shows the frictional force per unit logarithmic radius width A log a = 1 relative to the total 
frictional force, \j X B\/c per unit volume. Each panel is labeled with the density «h in cm" 3 , and 
each line is labeled with the grain charge q g . The solid lines are when the frictional force is parallel 
to the magnetic force, and the dotted lines are when it is anti-parallel to the magnetic force, apart 
from the component perpendicular to j X B. The dot-dashed line in the panel of nn = 10 10 cm -3 
shows the sum of the frictional forces exerted by grains of charge — e and e. 
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Fig. 8.— Same as Fig. 7 but for B = 0.15, 



-34- 




Fig. 9. — Electric conductivities as functions of the mean density of the cloud, nn, for the same 
case as in Fig. 4. The thin lines are for conductivity cry along magnetic field lines, and the thick 
lines are for Pedersen conductivity ap for the case of B = B cr . Contributions of electrons, ions, 
and grains are shown by dotted, long-dashed, and short-dashed lines, respectively, and the totals 
are shown by solid lines. The line labeled cth represents absolute values of Hall conductivity for 
the case of B = B CT ; the dot-dashed line is where cjh < 0, and the solid line with an > 0. 



